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ABSTRACT 


A  general  computer  program  for  the  calculation  of  the  chem¬ 
ical  equilibrium  composition  of  a  gaseous  system  was  written  for 
the  IBM  7090/94.  The  program  is  based  on  the  minimization  of 
the  Gibbs  free  energy  of  the  system;  the  resulting  non-linear 
equations  are  solved  by  a  modified  Newton-Raphson  iteration 
scheme.  The  three  options  presently  available  for  the  two 
intensive  variables  necessary  for  the  calculation  of  the  equi¬ 
librium  composition  are  pressure-temperature,  pressure-enthalpy 
and  pressure-entropy,  although  other  options  for  any  two  inten¬ 
sive  variables  may  be  readily  added. 

A  "program  generator  that  produces  the  source  program 
cards  in  FORTRAN  IV  of  the  equilibrium  p:.o-.,rum  for  a  specific 
chemistry  system  was  also  written.  in  the  process  of  gener¬ 
ating  the  source  program,  the  generator  extracts  from  a  library 
magnetic  tape  the  necessary  curve  fits  for  the  species  enthalpy 
and  entropy  as  functions  of  temperature  for  the  specified 
chemistry  system.  This  resulting  source  program  may  be  utilized 
as  a  subroutine  to  serve  the  needs  of  each  particular  application. 
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NOMENCLATURE 


the  number  of  atoms  of  element  j  in  one  molecule 
of  species  i 

the  left-hand  side  of  the  element  conservation 

matrix  [Eq.  (6)1  consisting  of  the  matrix 

elements  a. . 

:ij 

partitioned  A  matrix  consisting  of  NT  rows  and 
NP  columns 


partitioned  A  matrix  consisting  of  (NSP-NP)  rows 
NP  columns 


element  conservation  constants  from  Eq.  (5) 
+  tn  P  [Eq.  (9) ] 


hi 


T 


s . 
1 


the  Jacobian  (3P./SL  X  )  defined  bv  Eq.  (31) 

D  n  '*■' 

residual  vector  defined  by  Eq.  (30) 


Gibbs  free  energy,  cal 

fhi  1  -j 

-  s^J  +  -tn  P  +  tnX^j  [see  Eq.  (16)] 


/I  fh 
IR 


the  row  vector  consisting  of  the  g  [see  Eq,  (18)] 

j. 

the  partitioned  row  vector  G  consisting  of  NP 
vector  elements 

the  partitioned  row  vector  G  consisting  of  (NSE--NT 
vector  element  s 

mixture  enthalpy,  ft2/sec2 

mixture  enthalpy,  cal 

species  enthalpy,  cal/mole 


V 


h. 

1 

NP 

NSP 

P 

R 

A 

S 

S 

s 

s . 

1 

s’. 

l 

T 

X. 
l 

Y. 

l 

Y 

Greek 

$ 


n . 
J 


species  enthalpy,  ft3/sec2 

number  of  elements  (including  charge,  if  applicable) 
of  the  chemistry  system  considered. 

total  number  of  species  in  the  system 

pressure,  atmospheres 

universal  gas  constant  1.98726  cal/mole  °K  or 
89506  fta/sec3“°K  lb/lb  mole 

mixture  entropy  lb  mole/lb 

mixture  entropy,  cal/* K 

mixture  entropy,  ft2/sec2  °K 

species  entropy  at  standard  state  (p  =  1  atm) , 
cal/mole  °K 

species  entropy  at  standard  state  (p  =  1  atm) , 
ft2 /sec3  ° K 

temperature,  °K 

mole  fraction  of  species  i 


number  of  moles  of  species  i 
the  total  moles  in  the  system 


NSPTN'IP 

~  +  )  )  Y.a. .  -  b .  ft .  [see  Eg.  (8)] 

RT  L  L  l  11  3  3  ^  J 

j=l  Li=l  J 

Lagrange  multipliers  [see  Eq.  (8)] 

matrix  elements  resulting  from  the  matrix  multi¬ 
plication  AaAi1  [see  Eq.  (26)] 


the  column  vector  co"sisting  of  the  modified 
Newton-Raphson  correction  terms  [Eqs.  (39)] 


Kronecker  delta 
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I.  INTRODUCTION 

The  problem  of  computing  the  equilibrium  composition  of  a 
gas  has  been  the  subject  of  many  papers.  An  extensive  search 
of  the  literature  has  revealed  that  two  approaches  have  been 
used: 

1.  Utilization  of  equilibrium  constants  associated 
with  the  pertinent  chemical  reactions. 

2.  Minimization  of  the  Gibbs  Free  Energy  of  the  system. 

The  equilibrium  constant  approach  is  the  most  common  for 

small  systems,  where  the  system  cf  equations  is  comprised  of 
the  element  conservation  equations  and  the  laws  of  mass  action. 
Where  onlv  a  few  species  are  involved,  a  set  of  reactions  can 
be  easily  formulated,  for  which  curve  fits  of  the  corresponding 
equilibrium  constants  as  functions  of  temperature  can  be 
obtained.  In  recent  years  computer  programs  were  written  to 
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formulate  the  reactions  for  large  chemistry  systems.  Curve 
fits  of  the  species  Gibbs  Free  Energy  as  a  function  of  temper¬ 
ature  are  obtained,  and  the . equi librium  constant  is  computed 
from  the  change  in  the  Gibbs  Free  Energy  across  each  reaction 
at  a  reference  state. 

For  the  second  approach,  the  Gibbs  Free  Energy  of  the 
entire  system  is  minimized  subject  to  the  constraint  of  charge  and 
element  conservation.  A  variety  of  techniques  have  been  devel¬ 
oped  to  solve  the  resultant  system  of  non-linear  equations?  the 
numerical  technique  described  in  this  report  is  a  modified 
Newton  Raphson  iterati  >n  scheme. 

The  objectives  of  the  effort  described  herein  are  two¬ 
fold: 

1.  Create  an  equilibrium  program  which  may  be  used  as 
a  subroutine  to  generate  the  chemical  equilibrium 
composition  and  f'ssired  thermodynamic  properties 
for  any  specified  gaseous  chemistry  system. 

2.  Write  a  "program  generator"  which  automatically  pro¬ 
duces  the  &  jurce  cards  for  this  equilibrium  pro¬ 
gram  according  to  the  chosen  species,  which  are  input 
for  the  program  generator.  This  program  generator 


3 


extracts  the  necessary  thermodynamic  data  from  a 
library  magnetic  tape. 

Thus,  a  programming  system  has  been  created  which  provides 
the  user  with  an  analytical  tool  which  may  be  easily  linked  to 
computer  programs  designed  to  solve  a  wide  variety  of  fluid 
dynamics  problems.  The  equilibrium  program  generator  has  pro¬ 
duced  equilibrium  programs  which  have  been  successfully  used 
for  a  variety  of  applications. 

A  typical  application  of  such  a  generated  program  was  for 
the  computation  of  the  equilibrium  composition  of  an  Air-Teflon 
mixture  through  a  boundary  layer  adjacent  to  an  ablating  sur- 
ftee.  A  detailed  description  of  a  sample  application  is  shown 
in  Section  VIII. 
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II.  ANALYSIS 

The  Gibbs  free  energy  of  a  mixture  is  defined  as 

G  =  H  -  TS  (1) 

where  H  is  the  mixture  enthalpy  in  calories 

S  is  the  mixture  entropy  ...*  calories/0 K 
T  is  the  temperature  in  °K  . 

The  enthalpy  of  a  mixture  of  thermally  perfect  gas  can 
be  computed  from  the  species  enthalpy  for  a  given  temperature 
using  the  following  equation 


H=?Yihi(T)  (2) 

where  is  the  number  of  moles  of  specie  i,  h^(T)  Ls  the 
enthalpy  of  specie  i  in  calories  per  mole. 

The  entropy  of  a  mixture  of  thermally  perfect  gases  can 
be  computed  from  the  species  entropy  at  a  given  temperature 
and  pressure  using  the  following  equation 

S  -  S  =  L  Y . [  s  .  (T)  -  R  In  X.  -  R  In  Pi  (3) 

Oil  i 

where  P  is  the  mixture  pressure  in  atmospheres 

R  is  the  universal  gas  constant  in.  cal. /mole  °K 
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s^(T)  is  the  entropy  of  specie  i,  cal. /mole  ®K 

X.  is  the  mole  fraction  of  species  i 
1 

Sq  is  the  reference  entropy  taken  as  zero  at 
T  =  0  ®K,  p  =  1  atm. 


The  second  term  on  the  right-hand  side  of  Eq.  (3)  is  the 
contribution  to  the  entropy  due  to  the  diffusion  of  the 
species. 

The  third  term  is  the  contribution  to  the  entropy  for  an 
isothermal  process  where  the  pressure  is  chanced  from  one 
atmosphere  to  a  pressure  P. 

Substituting  Eqs.  (2)  and  (3)  into  Eq.  (1), 


G 

RT 


i-1 


s . 1  +  In  X .  +  In  P 
l  J  l 


(4) 


where  NSP  is  the  total  number  of  species  in  the  system. 

The  determination  of  the  equilibrium  composition  is 
equivalent  to  finding  the  set  of  X^'s  which  minimizes  (4), 
subject  to  the  constraint  of  the  element  conservation  equations 


NSP 

V  Y.  a . 

11 


i=l 


ID 


NP 


j~l»  2( « • . 


(5) 


6 


where  a. .  is  tne  number  of  atoms  of 
ID 

of  specie  i,  is  the  total  number 
NP  is  the  number  of  elements  in  the 
Written  in  matrix  form  Eq.  (5) 


element  j  in  one  molecule 
of  atoms  of  element  j ,  and 
system, 
becomes 


,Ya'**iYNSP] 

al  1 

ai  a 

•  •  • 

al  ,NP 

aa  i 

• 

aa  a 

•  o  • 

a  a ,  NP 

• 

aNSP#-j 

aNSP,e 

•  •  • 

aNSP,] 

(b^ba  o  c  •kjjp) 


(6) 


One  additional  equation  which  must  be  satisfied  by  the 
X . 1 s  is  that 


NSP 


i=l 


Eq. 


Eq  (4)  is  minimized  subject  to  the  constraint  given  by 
(5)  using  Lagrange  multipliers.  Let 


NP 


$>  =  --  +  \ 

RT  L 

3=1 


'NSP 

)  Y .  a  .  .  -  b 
L  i  n 


Li=l 


V 


(8) 


7 


(12) 
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Note  that: 


a  /  Y  \  Y  -  y 

Y£  (tn  ”y  )  =  Y 


for  l  =  k 


(13) 


V 


for  l  ^  k 


From  which  it  can  be  seen  that 


Y.  ) 

£.Y. 

C,  +  -tn 

i  l  dY. 

k 

k  1 

Y  i 

(14) 


If  Eq.  (14)  is  substituted  into  Eq.  (ll)and  this  result 
is  substituted  into  the  expression  for  the  partial  deriva¬ 
tive  of  [as  defined  by  Eq.  (8)]  with  respect  to  Y^f  then 


dY. 


NP 

C .  +  InX  .  +  )  a.  .  tt  . 

l  i  L  id  3 

j=l 


(15) 


It  can  be  shown  that  the  solution  —  =  0  represents  the 

dYi 

minimum  of  the  Gibbs  free  energy,  which  corresponds  to  the 
unique  solution  where  the  system  is  in  chemical  equilibrium 
(see  Ref.  l) . 

There  are  NP  element  conservation  equations  given  by  (5) 
ana  NSP  equations  obtained  through  the  minimization  process 


given  by  (15) .  The  unknowns  in  this  system  of  equations  are 
NSP  species  mole  fractions  and  NP  Lagrange  multipliers. 
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The  system  of  (NSP  +  NP)  equations  with  the  same  number 
of  unknowns  will  be  reduced  to  a  system  of  NSP  equations  with 
NSP  unknowns. 

Let 

h. 

g.  =  -  C.  -  tn  X.  =  -  -  s.j+  In  P  +  In  X.} 

(lo) 

then  Eq.  (15)  can  be  written  as 

NP 

y  a.. 77.  =  g.  i=l,  2,  . . .  ,NSP  (17) 

L  ij  j 

j=l 


for  the  condition  of  the  minimum  Gibbs  free  energy  (i.e. 


=  0  . 


If  Eq.  (17)  is  written  in  matrix  form  the  result  is 


(18) 
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or 


A  n  =  G 


(19) 


The  set  of  equations  represented  by  Eq.  (18)  can  be  written 
as  one  set  of  NP  equations  and  another  set  of  (NSP-NP)  equa¬ 
tions.  If  Eq.  (19)  is  partioned  accordingly  the  result  is 


or 


NP  M 

i , 

NSP-NP 


NP 


NSP-NP 

(20) 


(21) 


and 


A3n  *  Ga  (22) 

Solving  Eq.  (21)  for  T1  =  A:  1  Gj,  and  substituting  into  Eq.  (22), 
one  obtains: 


A3  a^  g^  — 


(23) 


If  the  expression  for  g^  given  by  (16)  is  substituted 
into  (23)  the  result  is 


In  X 

n 


NP 

l  {^n-NP, j 

j=l 


+  In  X . )\  -  C  (24) 

3  J  n 

n  =  NP  +  1, NP  +  2, . . . 


,  NSP 


11 


C.  =  - 
1 

g ,  -  £n  X . 

_  1 

R 

f  h.  1 

[T  "  si_ 

+  in  P 

111 

$1  a 

$1  3 

*•*  ^ l , NP 

$3  1 

$3  3 

$3  3 

***  ^a,NP 

[  Aa  Ax  ]  - 


$ 


$ 


NSP-NP,  l  NSP-NP, a  NSP-NP, 3 


$ 


NSP-NP, NP 


(26) 


Eqs.  (5)  and  (24)  correspond  to  NSP  equations  in  the  NSP 
unknowns  (mole  fractions).  In  addition  Eq.  (7)  can  be  used 
as  a  check  equation.  However,  since  the  total  number  of  moles 
in  the  system,  may  be  of  interest,  this  parameter  is  added  to 
the  set  of  unknowns  and  Eq.  (7)  is  added  to  the  system  of 
equations . 

This  is  accomplished  by  writing  Eq.  (5)  in  the  following 

form 


NSP 


)'  X.  a.  . 

L  i  ID 

i=l 


NP 


Y 


D  1 ,  .  .  .  , 


(27) 


12 


where  Y  is  the  total  number  of  moles  in  the  system,  i.e. 


NSP 

Y  =  V  Y. 

Li  1 

i=l 


X. 

1 


(23) 


(29) 


The  non-linear  system  of  equations  given  by  (7) ,  (24) 
and  (27)  in  the  unknowns  Y  and  X^  (i=l, 2, . . .NSP)  is  solved  by 
using  a  modified  Newton-Raphson  technique. 
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III.  METHOD  OF  SOLUTION 


A  residual 

vector, 

represented  by  F_.  , 

is  obtained  from 

Eqs.  (7)  and  (27) 

as  follows: 

NSP 

F>  “  l  xi 

i  =  l 

-  1 

F 

NSP 

=  )  X  •  a  .  . 

D  L  ii] 

i=l 

. 

Y 

(30) 

j=2, 3, . . 

o ,  NP 

The  Jacobian  to  be  used  in  the  Newton 

Raphson  iterative  pro- 

cedure  for  NP  dimensions, 

DF  =  (3FV3 

In  X  ) 

for  (1  ^ 

j  *  NP)  , 

(1  £  l  *  NP)  is  formed  as 

follows 

dFj 

3FX 

3FX 

3  In  Xx 

3  In  X3 

•  •  • 

3  -tn  X 

NP 

3Fa 

3Fa 

3Fa 

3  In  Xx 

3  -tn  Xa 

•  •  • 

3  XNP 

DF  = 

• 

• 

• 

• 

(31) 

• 

3F 

NP 

• 

3F 

NP 

• 

3fnp 

a  In  Xx 

3  -tn  Xa 

•  •  • 

»  XH£_ 

14 


ious  terms 


The  following  is  required  for  the  evaluation  of  the  var- 


5  4 


3  In  X . 

•J  _  c 

3  £n  X^  j  t 


sxj  _ 

3  In  X. 


/ 

1  ^  j  s.  NP 

\  /  (32) 

1  S  (  £  NP 

y 


where  6 . ,  is  the  Kronecker  delta. 
jl 


1  for  j  =  l 


) 


V"  < 


0  for  j  ^  t  J 


(33) 


From  Eq.  (24)  for  (NP  +  1  £  n  £  NSP) ,  (1  *  l  s  NP) 


3  In  X 


n 


NP 


3  In  X,  3  In  X. 

-L 


I  [(Ck  +  V*n-HP.k]-Ca 


k-1 


and 


n-NP,  l 

(34) 


3X 


n 


=  x  8 

d  (n  X.  n^n-NP, l 


since  Ch  -  Ch(P,T). 
From  Eq.  (27) 


1  _  i 

Y  ~  bi 


NSP 


)  X  .  a  . 

L  ii- 


i  =  l 


(35) 
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then 


&  In  Xl 


Xlali+ 


NSP 


V  X.a  /3  f 

L  j  ni  n-NP,£ 

n=NP+l 


(36) 


It  should  be  noted  that  Sqs.  (35)  and  (36)  are  valid 
only  for  1^  /  0. 

The  partial  derivatives  that  are  required  for  DF  as 
defined  in  Eq.  (31)  can  now  be  evaluated  as  follows: 


3F 


NSP 


B  In 


hrmXt*  1  xA-np, 


n=NP+: 


NSP 


/  \ 
1 


a  -tn  X. 


>  [X.a  .0  XTn  - 
L  1  ni  n-NP,l 


X.a.,  + 

l  Ij  l,  j  nj 

n=NP+l 


b. 


(37) 


j  b  In  Xt) 


j-2 , 3 , . . . ,NP 
, 2 , .  .  . ,NP 


If  the  residual  vector  given  by  Eq.  (30)  is  denoted  by  the 
column  vector 


Fi 

F3 


(38) 
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and  the  set  of  correction  terms  denoted  by  the  column  vector 


A  tn  Xx 
A  -In  Xa 


9 


A  Ia 


(39) 


and  DF  by  Eq.  (31) ,  then  in  accordance  with  the  Newton- 
Raphson  procedure  these  correction  terms  are  found  by  solving 
the  following  set  of  simultaneous  linear  equations 


(DF)  *  (A)  =*  F 


(40) 


The  next  set  of  iterated  values  of  the  mole  fractions 
of  the  "prime"  species  is  given  by  the  expression 


{ln  xjW  =  {ln  Vn  ■ X(A  Vn  (41) 

j=l, . . . ,NP 

where  N  denotes  the  iteration  number  and  X  is  a  relaxation 
factor  (0  <  X  s  1) • 

The  relaxation  factor  is  chosen  by  the  program  so  that 
it  is  as  large  as  possible  (but  never  larger  than  unity)  and 
such  that  ma:c  {|f.||  decreases  from  iteration  to  iteration. 
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Initially  X  assumes  a  value  of  unity,  but  is  subsequently 
reduced,  if  necessary,  until  max  j j F j  jj  is  smaller  than  max 


{!■ 


'N 

The  relaxation  factor  X  has  been  introduced  to 


N-l 


keep  the  iterations  in  the  neighborhood  of  the  solution  when¬ 
ever  possible.  If  the  solution  is  converging  and  X  <  1,  then 
X  is  increased  subject  to  the  above  constraint  to  accelerate 
convergence. 

The  mole  fractions  for  (NP  +  1  £  n  £  NSP)  are  obtained 
from  Eq.  (24) ,  and  the  total  moles  in  the  system  from  (35) . 

The  criteria  ior  convergence  of  the  iteration  procedure 


are  as  follows: 


(1) 


F.  j  <  10"6,  j=l, 2 , . . . , NP 


where  the  F_.  are  defined  by  Eq.  (30)  ,  and 


(2) 


(to  x.)N  -  (to  x.)N-1 

(t“  Vn 


'  10 


-4 


i=l, 2,3, . . . ,NSP 


where  N  is  the  iteration  counter. 


The  first  test  ensures  that  the  major  species  have  con¬ 
verged,  and  the  second  test  applies  to  all  the  species  in  the 
system  including  the  trace  species,. 
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I V .  SOME  REMARKS  ABOUT  CONVERGENCE 

Mathematical  analysis  (s-e  Ref.  1)  has  shown  that  the 
Gibbs  free  energy  function  is  strictly  convex;  hence,  setting 
the  derivatives  of  Gibbs  free  energy  to  zero  results  in  a 
unique  solution.  However,  since  the  system  of  equations  can¬ 
not  be  solved  in  closed  form,  the  question  arises  as  to  under 
what  conditions  will  the  program  encounter  convergence 
problems. 

It  is  suspected  that  the  majority  of  convergence  problems 
will  be  cittributable  to  either 


the  estimates  of  the  mole  fractions  Xy  j=l, 2, . . . ,NP, 
are  too  far  from  the  solution, 


(2)  the  species  have  not  been  properly  ordered. 

The  user  must  consider  the  ordering  of  species  when  using 
the  "program  generator."  The  first  NP  species  specified  must 
be  linearly  independent.  (Note  that  charge  conservation  is 
handled  exactly  like  element  conservation.)  For  example,  if 
a  seven  species  air  (Na ,  N,  03 ,  0,  NO,  N0+,  e“)  chemistry 
system  is  considered,  NP  =  3  (the  three  "elements"  are  N,  0, 
e~) ,  NSP  =  7  (there  are  7  species)  and  the  ordering  of  the 
first  three  species  could  be  Ns ,  0a ,  e”  or  N0+,  NO,  0a ,  as  well 


as  a  number  of  other  combinations.  However,  the  /-hoice  of  Na , 


N,  e”  as  the  first  three  are  incorrect,  since  N3  and  N  are  not 
linearly  independent. 

The  first  NP  species  must  be  chosen  such  that  each 
element  (and  charge,  if  applicable)  be  included  among  these 
species  This  is  a  necessary  condition  for  the  system  to  be 
linearly  independent;  however,  :  t  is  not  sufficient  to  guar¬ 
antee  linear  independence.  Consider  again  the  seven  species 
air  case,  and  let  the  first  NP  species  be  e”,  NO  and  N0+. 

If  one  forms  the  Ax  matrix  as  described  in  (20) ,  it  is 
apparent  that  this  system  is  linearly  dependent,  since  sub¬ 
tracting  the  second  row  from  the  first  results  in  the  third 
row.  However,  the  rule  described  above  will  work  in  the 
large  majority  of  chemistry  systems. 

Furthermore  the  program  will  work  most  efficiently  if 
the  first  NP  species,  besides  being  linearly  independent,  are 
chosen  so  that  their  mole  fractions  are  as  large  as  possible. 
The  user  should  also  note  that  the  species  cannot  be  ordered 
such  that  electrons  is  the  first  species,  since  Y  (total 
moles)  becomes  singular  as  computed  from  Eq.  (35) . 
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V .  OTHER  THERMODYNAMIC  OPTIONS 

As  described  in  the  previous  sections  the  program  computes 
the  equilibrium  composition  of  a  gaseous  mixture  for  specified 
values  of  pressure  and  temperature.  Other  input  options 
available  are  pressure-enthalpy  and  pressure-entropy.  For 
these  options,  the  program  performs  a  series  of  pressure  temp¬ 
erature  calculations  in  an  iterative  manner  until  the  specified 
values  of  enthalpy  or  entropy  are  satisfied  within  prescribed 
tolerances.  Estimates  of  the  corresponding  values  of  tempera¬ 
ture  and  species  mole  fractions  must  be  transmitted  to  the 
equilibrium  program  through  the  calling  sequence. 

Utilizing  the  specified  value  of  pressure,  and  the  esti¬ 
mated  values  of  temperature  and  the  first  NP  species  mole 
fractions,  the  program  computes  the  corresponding  equilibrium 

values  of  the  species  mole  fractions  (X^).  ,  ar^  of  enthalpy, 

—  /\ 

Hi  (or  entropy  Sl  ) .  ?  second  guess  of  the  temperature  is 

obtained  (for  enthalpy) : 

T2  =  Ti  i  0.05  Ti  for  Hx  £  H 


~  2  /  2 
where  H  is  the  specified  mixture  enthalpy  per  unit  mass,  ft  /sec*" 


A  S  lb  mole 

S  is  the  entropy,  —  ,  - rr - 

R  Id 


S  is  the  entropy  per  unit  mass, 
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The  program  then  calculates  Ha  corresponding  to  the 
specified  pressure  and  Ta .  The  third  estimate  of  temperature 
is  calculated  by  the  application  of  the  "regula  falsi"  method 
(or  method  of  chords),  as  follows: 


(43) 


All  subsequent  estimates  of  T  are  calculated  by  expres- 
sing  f^  as  a  function  of  T^  by  means  of  a  LaGrange  second- 
order  equation  and  evaluating  this  parabola  at  f  =  0,  to  obtain  the 

next  estimate  of  T^+^,  where  k  denotes  the  enthalpy  iteration  counter. 

H v  -  H  _6 

Convergence  is  assumed-when  -  £  10 

IT 

This  procedure,  of  course,  applies  also  to  the  pressure- 
entropy  option.  The  program  logic  may  be  readily  extended 
to  include  input  options  for  any  two  intensive  properties. 
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VI.  INSTRUCTIONS  FOR  USING  THE  GENERATED  EQUILIBRIUM  PROGRAM 


Initializing  Call  Statement 


CALL  GENL7X(X,,  PROP,  PR  ' ,  T,  P, KOPT, BEE) 


This  CALL  statement  initializes  the  addresses  of  the  formal 
parameters  in  the  calling  sequence  and  computes  the  element 
conservation  constants  described  by  Eq.  (27) ,  utilizing 
the  values  of  the  mole  fractions  stored  in  the  X  array,  upon 
entry.  The  value  of  KOPT  must  be  a  positive  integer.  These 
values  of  .he  mole  fractions  must  truly  represent  the  gaseous 
system  at  some  reference  state;  they  cannot  be  approximations 
(or  guesses)  of  the  X^. 

The  user  also  has  the  option  of  directly  storing  the 
element  mole  fractions  in  the  BEE  array  rather  than  specifying 
the  X^'s.  The  option  of  specifying  the  values  directly  in 
the  BEE  array  is  particularly  useful  for  applications  where 
the  element  mole  fractions  are  not  constant,  such  as  the  case 
of  a  boundary  layer  adjacent  to  an  ablating  body.  When  the  con¬ 
tents  of  the  BEE  array  are  specified  directly,  KOPT  must  be 
a  negative  integer  upon  entry,  /'ll  subsequent  calls  (see 
below)  must  then  specify  KOPT  <  0,  and  the  BEE  array  must  con¬ 


tain  the  current  values  of  b_.  upon  entry. 
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The  formal  parameters  in  the  calling  sequence  are  described 
below: 

X  is  a  one-dimensional  array  which  contains  the  values  of 
the  mole  fractions-  The  size  of  the  array  is  equal  to 
one  more  than  the  number  of  species  in  the  system  (NSP+1). 

o 

T  is  the  temperature  m  K. 

P  is  the  pressure  in  atmospheres. 

KOPT  is  a  code  word  which  specifies  the  input  option 
( see  below) . 

BEE  is  a  one-dimensional  array,  whose  size  is  equal  to 
the  number  of  elements,  NP,  which  contains  the  element 
conservation  constants  b.  described  by  Eq.  (27). 

J 

PROP  and  PR 2  are  dummy  variables  whose  contents  are 
identified  by  the  code  word,  KOPT  (see  below). 


General  Call  Statment 

CALL  GENR7X 

This  CALL  statement  may  be  executed  repeatedly,  once  the 
initializing  CALL  statement  has  been  executed.  Note  that  the  equi¬ 
librium  subroutine  operates  upon  those  formal  parameters  which  appear 
as  arguments  in  the  initializing  CALL  statement.  Hence  these 
arguments  must  contain  the  correct  values  of  the  input  param¬ 
eters  when  this  general  call  statement  is  executed.  It  is 
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suggested  that  these  formal  parameters  (X,  PROP,  PR2.  T,  P, 

XOPT  and  BEE)  be  included  in  COMMON  storage  so  that  transmis¬ 
sion  of  their  contents  among  all  routines  may  be  satisfactorily 
achieved. 


Pressure-Temperature  (P-T)  Option i  KOPT  =  ±  1 

"uesses  of  the  first  NP  species  mole  fractions  (X ^  , 
j=l, 2 , . . . , NP)  must  be  specified  in  the  X  array  upon  entry. 


Pressure-Enthalpy  (P-H)  Option:  KOPT  =  ±  2 


In  addition  to  these  guesses  of  X_.  ,  j=l ,  2 ,  .  .  . ,  NP,  a  guess 
of  temperature  must  be  specified  in  the  T  parameter  upon  entry. 


,  A , 

Pressure-Entropy  (P-S)  Option:  KOPT  -•  ±  3 

Same  inputs  as  for  the  P-H  option,  above. 

KOPT  =  ±  1  for  pressure-temperature  input  option 
±  2  for  pressure-enthalpy  input  option 
±  3  for  pressure-entropy  input  option. 

If  KOPT  >  0  the  contents  of  the  BEE  array  (values 


of  b_.) 


are  not  altered. 

If  KOPT  <  0  the  user  must  specify  the  values  of  b .  in  the 


array  BEE  upon  entry. 
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1,^  contents  of  PROP  and  PR2  are  described  in  the  table 
below.  For  KOPT  =  ±  1,  the  contents  of  PROP  and  PR2  are  output 
properties;  for  KOPT  =  ±  2  and  ±  3,  the  content  of  PROP  is  a 
specified  input  and  the  content  of  PR2  is  an  output  property. 


KOPT 

PROP 

PR2 

± 

1 

enthalpy, 

f t2/sec2 

entropy, 

lb  moles/lb 

± 

2 

enthalpy, 

ft2/.«ec2 

entropy, 

lb  moles/lb 

4- 

3 

entropy, 

lb  moles/lb 

enthalpy, 

f t2/sec2 

There  are  four  situations  under  which  the  equilibrium  pro¬ 
gram  will  call  an  error  subroutine,  and  they  are: 

1.  The  temperature  at  which  the  species  enchalpy  and 
species  entropy  are  to  be  evaluated  from  the  curve  fits,  as 
shown  in  Eqs.  (44) ,  is  higher  than  the  upper  bound  of  tempera¬ 
ture  range  over  which  the  fit  is  valid  (no  check  is  made  for 
the  lower  bound  of  the  temperature). 

2.  If  more  than  ten  iterations  are  required  to  satisfy 

the  prescribed  value  enthalpy  (or  entropy)  within  a  relative 
—  0  —  /\ 

error  of  10  ,  when  the  P-H  (or  P-S)  options  are  used 

3.  If  more  than  30  iterations  are  required  to  converge 
to  a  solution  for  a  P-T  problem. 
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4,  If  between  two  successive  iterations  for  a  P-T  problem, 
max  |  Fj  j j  as  defined  by  Eq.  (30)  does  not  decrease,  the  value 
of  the  relaxation  factor  >,  used  in  Eq.  (41)  will  be  reduced; 
this  r  duction  of  X  will  be  repeated  until  max  |JFj  |j  decreases 
from  iterati.cn  to  iteration.  However,  the  program  calls  the 
error  subrout' ne  if  Xbas  been  reduced  for  10  consecutive  P-T 
iterations . 

For  the  first  three  situations  it  will  probably  be 
desirable  to  let  the  program  continue  after  writing  an  error 
message.  If  tne  last  situation  occurs,  an  error  message  should 
be  output  after  which  the  execution  should  be  terminated.  The 
run  should  be  resubmitted,  either  with  "better"  guesses  for 
the  mole  fractions,  or  possibly  the  species  may  have  to  be 
re-ordered. 

A  listing  of  the  error  subroutine  which  handles  these 
situations  as  stated  above  appears  in  Appendix  c  • 

Included  in  the  equilibrium  package  is  a  subroutine  that 
will  solve  a  system  of  N  simultaneous  equations.  The  method 
employed  is  pivotal  condensation,  where  the  largest  pivotal 
element  available  is  utilized. 

The  user  may  CALL  this  subroutine  in  another  part  of 
hi.s  program  by  the  following  statement: 


27 


CA_iL  CLEM  (Al,  DLNX,  F,  N,  N+l,  AT) 

where  the  formal  parameters  are  defined  as  follows: 

Al  is  the  matrix  of  coefficient 
DLNX  is  the  solution  of  vector 
F  is  the  forcing  vector 

N  is  the  number  of  equations  to  be  solved 
AT  is  a  working  array  used  by  the  subroutine,  but 
must  be  dimensioned  by  the  calling  program  to  be 
of  size  (N,  N+i) . 


Streamline  Calculation  Using  Equilibrium  P-S  Option 
To  perro^m  calculations  along  streamlines  using  the  P-S 


option  of  ^he  equilibrium  chemistry  program,  it  is  first 
necessary  to  execute  either  the  P-T  option  or  the  P-H  option 


at  the  initial  point  on  the  streamline  in  order  to  extract 
the  value  of  entropy  (S).  All  subsequent  calculations  along 

the  streamline  may  then  be  executed  using  the  P-S  option  where 

A 

S  is  the  value  extracted  at  the  initial  point.  This  value  of 

A 

S  is  conserved  along  each  streamline.  To  calculate  the  dimen- 

A 

sionless  value  of  entropy  it  is  necessary  to  multiply  S  by  the 


mixture  molecular  weight. 
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VII .  DESCRIPTION  OF  THE  EQUILIBRIUM  PROGRAM  GENERATOR 

For  each  chemistry  system  used  the  user  must  first  gen¬ 
erate  the  FORTRAN  IV  source  cards  for  the  specific  chemistry 
system  requested.  The  generator  program  punches  the  source 
cards  of  the  equilibrium  subroutine,  including  all  the 
required  internal  data,  and  the  elements  of  the  A  matrix 
defined  by  Eq.  (6) .  The  curve  fits  are  of  the  form: 


c 


Pi 


a  .  +  b .  T 
1  i 


_  T  b . T2 

h,  =  h .  +  f  C  dT  =  a .  T  +  +d.  (44) 

l  l  J  p .  \  2  l 

o  o  i 


a.  £n  T  +  b .  T  +  c . 
i  ri 


[T>T 

i  ° 


i 


where 


is  the  species  specific  heat  at  constant  pressure, 
ft2/sec2  °K 


h.  is  the  species  enthalpy,  ft2/sec2  (note  that 

species  sensible  enthalpy  is  taken  to  be  zero  at  0  K) 

h.  is  the  species  heat  of  formation  at  0  °K,  ft2/sec3 
i 

o 

s.  is  the  species  entrcpy,  ft2/sec2  °K  (note  "s .  is 
taken  co  be  zero  at  0  °K,  p  =  1  atmosphere) 


T 


o . 
i 


is  the  lower  bound  of  temperature  for  the  curve  fit 

o . , 

or  c  ,  K. 
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The  data  for  Cp^  have  been  fit  (in  the  3  sast  square 
sense)  as  a  function  of  temperature,  with  first  order  equa¬ 
tions;  the  curve  fits  are  subdivided  into  several  temperature 


ranges  The  only  constraint  on  these  fits  for  the  equilibrium 
program  is  that  the  temperature  ranges  for  each  fit  (rip  ,  h^, 
s.)  of  a  given  species  coincide.  The  curve  fits  now  con¬ 
tained  on  the  library  tape  used  by  the  program  generator 


satisfy  this  constraint.  Presently  a  maximum  of  6  tempera¬ 


ture  ranges  per  species  is  accommodated. 
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Input  Format  to  the  Program  Generator 


CARD  NO. 


2 


COLUMNS 


3-8 


45 

46 

1 

43-45 


3 


1 


4,5,6,...  6A12 


DATA  TO  BE  PUNCHED  ON  CARD 

Date  (month/day/yea r ) 

1 

1 

7 

(right  adjusted)  number  of  species 
in  system 
2 

The  species  symbols  consistent  with 
GASL  library  tape,  with  12  columns 
per  field  for  each  species.  The 
first  NP  (where  NP  is  the  total 
number  of  elements  and  charges,  if 
applicable)  species  must  be  linearly 
independent,  and  the  first  of  these 
species  may  not  be  electrons  since 
Eq.  (35)  would  become  singular. 


See  Appendix  D  for  sample  inputs 


Operating  Instructions  for  the  Program  Generator 

When  running  the  "program  generator"  on  a  7094  IBM  computer 
the  following  operating  instructions  ripply : 

1.  IBSYS  system  version  13  is  to  be  used. 

2.  Mount  a  specified  GASL  Chemistry  Library  tape  on 
unit  A5. 

3.  Scratch  tapes  are  to  be  mounted  on  units  B5  and  B6. 

4.  One  file  from  output  tape  unit  A6  is  to  be  punched 
on  FORTRAN  cards;  interpret  columns  1-60,  These 
cards  are  the  source  deck  of  the  equilibrium 
program. 


5. 


List  under  program  control  the  output  from  unit  Bl„ 
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Output  Obtained  from  the  Equilibrium  Program  Generator 

The  generator  will  output  (on  paper)  a  report  detailing  the 
coefficients  of  the  least-square  curve  fits  of  the  species 
enthalpy  and  species  entropy,  as  functions  of  temperature  as 
well  as  the  upper  and  lower  bounds  of  the  temperature  regions 
corresponding  to  each  curve  fit.  Since  the  equilibrium  pro¬ 
gram  does  not  check  for  the  lowest  bound  on  temperature  (i.e. 
it  will  automatically  extrapolate  to  a  temperature  below  the 
lowest  bound) ,  the  user  shpuld  check  this  output  if  he  expects 
the  program  to  encounter  temperatures  below  the  lowest  bound, 
which  for  most  species  presently  accommodated,  is  200°K. 

In  addition,  the  source  deck  of  the  equilibrium  sub¬ 
routine,  in  FORTRAN  IV,  for  the  specified  chemistry  system  is 
punched  as  output,  a  listing  of  which  appears  in  the  appendix. 
The  user  must  add  subroutine  ERROR,  as  described  in  Appendix  C. 


VIII . 


SAMPLE  APPLICATION 


A  typical  application  of  a  generated  equilibrium  program 
was  for  the  computation  of  the  equilibrium  composition  of  an 
Air-Teflon  mixture  through  a  boundary  layer  adjacent  to  an 
ablating  surface.  An  equilibrium  subroutine  comprised  of  37 
species  was  produced  by  the  "program  generator,"  A  "user's" 
program  was  written,  utilizing  this  subroutine,  which 
traversed  the  boundary  layer,  calculating  the  equilibrium 
composition  of  the  gas,  and  the  mixture  parameters,  such  as 
the  temperature,  corresponding  to  a  specified  variation  of 
enthalpy,  pressure,  and  element  composition.  The  element 
composition  varied  from  90%  teflon  to  air  ratio,  to  100%  air, 
and  the  temperature  varied  from  1000  °K  to  5500  °K.  The 
frontispiece  depicts  th-  variation  of  the  mole  fractions  with 
temperature  for  a  preliminary  chemical  system  consisting  of 
28  species  which  reflects  a  teflon  air  ratio  of  50%  by  volume. 
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APPENDIX  A 

FLOW  CHARTS 
OF  THE 


EQUI LIBRIUM  PROGRAM 


FLOW  CHART  FOR  SUBROUTINE  GENL7X 


ENTRY  GENR7X 


<0 


I 


Test 

KOPT 


s0 


Store  bee(i) 

In  A (NSPl / I ) 
1  ^  I  ^  NP 


KOPT=-KOPTj 


Do,P,T  problem  (see  separate  flow  chart 
starting  at  Q 


\  >10 

'-■”1 


I  CALL 

i  TTDPOD 


Return 


FLOW  CHART  FOR  P,T  PROBLEM 


Compute  Fj (1  ^  j  ^  N?)  eq.  30 

Store  Max  |f j  |  -*  F  (NP  +  1) 

Store  F,*  FOLD (I)  (1  s  j  s  NP  +  1) 

(B  .  . 


Compute  Jacobian  eq.  31,  Solve  for  At-nX^  eq.  40 


Compute  £  i  s  NSP)  eqs.  41,  24  and  store-*XLN(I) 

Compute  (1  ^  i  s  NSP),  Any  X^  >  1  set  tnXi  -  -.COi 


APPENDIX  B 

LISTING  OF  A 
GENERATED  PROGRAM 


ete) 


O'JOOO^OOO 


B-l 


SIBHC  GEN7X.  FUIIST 

SUBROUTINE  GEM  7XiX,  PROP.  PK2,  T  »  P,  KOP  T  » BEE) 

DIMENSION  AH  4,  4), XI  9 1 « BE  TA  (  4,  4),Ai  9,  4),E(  5),Zt  B), 

IXLM  8  )  »  C  I  8 )  ,XLAST I  81 ,ULNX<  4),PSlt  4,  4,  4), 

2DENOMI21 , TEMPI  4), TEMPI!  4), ATI  4,  5),FQIC(  5), 

3N  <  8), EMI  8 ) t CP  1  8,6,5) , BEE  I  4) 

SPECIES  CROER 

1  C2 

2  N2 

3  £- 

4  AR 

5  N 

6  NC 

7  )iC  + 

8  C 


DATA  fe  / 


I 

2.798125CE 

C3, 

3.1960308E 

03, 

1.632C951E 

C8, 

2.2416383E 

03, 

2 

6.3920616E 

C  3 , 

2.98387  ICE 

03, 

2.S838710E 

03, 

5. 596250GE 

03/ 

DATA  EM 

/ 

1 

3.20CC00CE 

Cl, 

2 • 80 160CCE 

01, 

5.4862000E- 

04, 

3.9944 OOGE 

01, 

2 

1.4CC80CCE 

Cl, 

3.C0C800GE 

01, 

J.CCGBOOOE 

01, 

1 . 60C0000E 

Cl/ 

DATA  A 

/ 

1 

2.UCGOOCCE 

cc, 

0. 

, 

0. 

, 

0. 

, 

2 

C. 

f 

1 .CQCGOOCE 

oc, 

l.CCOCOOOE 

CO, 

l.COOOOOCE 

00, 

3 

0. 

f 

0. 

, 

2  •  GOO  OOOCE 

CO, 

0. 

, 

4 

0. 

f 

1 . OOOOCCOE 

00, 

l.COOOOOOE 

ou, 

1 .  OOOOOOOE 

CO, 

5 

c. 

f 

0. 

• 

C. 

, 

0. 

, 

6 

I.OCCOOOCE 

cc, 

0. 

f 

0. 

C. 

, 

7 

-  I • CCCOGOCE 

CO, 

0. 

, 

0. 

, 

0. 

, 

b 

0. 

f 

0. 

, 

l.CCOCOOOE 

00, 

0. 

, 

9 

0. 

9 

0. 

, 

c. 

, 

G. 

/ 

DATA  BETA 

/ 

1 

-0. 

, 

5 . OOCOOOOE 

-01, 

5.CCCG000E- 

01  , 

5.  COCGOOUE- 

-ci, 

2 

5 • COGOOOOE 

-01, 

5 . OOCCCOCE 

-01, 

5.CCCOOOOE- 

01  , 

0. 

, 

3 

0. 

» 

0. 

- 1  • OOOOOOGE 

CO, 

0. 

, 

4 

G. 

f 

0. 

, 

0. 

, 

0. 

/ 

DA7AICP 

I  III 

,11,1=  1 

, 

76)/ 

I 

9.3716999E 

C3, 

1 . C  7  30000  E 

04, 

4 .C757999 £ 

C  8 , 

5. 5S25999E 

C3, 

2 

1.S922000E 

C4 , 

1 .01 7300CE 

04  , 

1 .CC15C00E 

04, 

1 • 61 BdOOCE 

C4 , 

3 

1. 1342000E 

04, 

1 .3349000E 

04, 

4.C7579S9E 

C8, 

1. 7240000E 

03, 

4 

1.18C200QE 

C4 , 

9 . 3989999t 

03, 

1  •  2823000E 

C4 , 

1.415400CE 

C4 , 

5 

1 • 25  75000E 

C4 , 

l  .3870GC0E 

04 , 

4.C757999E 

C8  , 

-1.284500UE 

03, 

6 

2.015G0CCE 

04, 

1 .245000CE 

04, 

6.8981000E 

03, 

1 . 3418000E 

C4 , 

7 

I.6IC800GE 

04, 

6 .2 1 57000E 

03, 

4.C757999E 

C8 , 

- 1 • 3CC900GE 

04, 

8 

-I.5187CCCE 

04, 

1  .2587CC0E 

04, 

1 .  1969000  E 

C3  , 

1.44870006 

03, 

9 

1.765200Gb 

04, 

2.628300CE 

C  3 , 

4.C757999E 

Cb , 

-2.488300CE 

04, 

1 

2.22410CCE 

C4 , 

1 .62250C0E 

04  , 

2.222 1G0GE 

03, 

-1.8841000E 

04, 

2 

2.C238GC0E 

C4 , 

1 .2720GC0E 

04, 

4.C757999E 

Cb  , 

-5.11 1900GE 

C4 , 

3 

5.565100CE 

C4  , 

2 .2964CC0E 

04, 

2.283 1000E 

04, 

7. 96 17GG0E 

03, 

4 

2.308500CE 

CO, 

1 .7330U00E 

CO, 

0. 

, 

2.7274GCGE 

-03, 

5 

5. 1C45000E 

-02, 

l .39C50CCE 

00, 

1 • 5  749CC0E 

CO, 

-4. 135300GE 

CO, 

6 

6.6608959E 

-Cl, 

2 •  77  76CCCE 

-01, 

0. 

, 

3. 5441CC0E 

-Cl, 

7 

* 

1.299500CE 

cc, 

2.681 100CE 

00  , 

1.7045C00E- 

01  , 

-6 • 72  38000E 

-02, 

8 

3.85890CCE 

-Cl, 

1 .37C10C0E 

-01, 

0. 

, 

6 • COG 100GE 

-01, 

9 

3. 174800CE 

-01, 

3.3365CCCE 

-01, 

1 • C 1690COE 

CO, 

2.27110C0E 

-Cl, 

1 

-<.C1760CCE 

-02, 

1 . C375C0CE 

CO, 

0. 

, 

1.468500CE 

00/ 

UATAICP 

(1,1 

,i),I  =  77 

,  152)/ 

1 

2.52600GCE 

cc, 

2 . 88 19CCCE 

-01, 

1.617 10GOE 

CO  , 

l.Cdb-'OOGE 

CO, 

2 

-l .48870CCE 

-Cl, 

1 .J7920CCE 

CO, 

0. 

f 

2.28740CCE 

cc, 

B-2 


3  9.983500GE 

-01, 

4.9663CCCE- 

-02, 

1.5316G00E 

CO, 

2.18 16000E 

CO, 

4-3. 13C30CCE 

-Cl, 

7.48460C0E 

-Cl, 

0. 

, 

3  •  9532000E 

00, 

5-2.39C1000E 

-Cl, 

-2.87290CCE 

-01, 

1.57760C0E- 

01, 

l. 15C8C00E 

CO, 

6  1 • 48  7500CE 

C4 , 

1.1777000E 

04, 

-1.9130000E 

C9  , 

9. 826 1999E 

03, 

7  2.6945000E 

C4 , 

1 .7289C0CE 

C4 , 

1 • 3522000E 

C4 , 

1. 72560006 

04 , 

8  2 . 8 742000E 

03, 

-5.238300CE 

C  3  , 

-1.913COOOE 

C9 , 

4. 1957000E 

C4 , 

9  5.62C30CCE 

C4 , 

2 • 1468000E 

04, 

-5.C189000E 

03, 

2.7862C0GE 

C4 , 

1-6 .2267U00E 

C3 , 

-8.99630COE 

C3  » 

-1 .913COOOE 

09, 

6. 7269000E 

C4 , 

2-1.C97500CE 

04, 

2 .64  C500CE 

03, 

4. 152C000E 

C4 , 

3.  28840CCE 

C4 , 

3-3.475000CE 

04, 

5.26C500CE 

04, 

-1.9130000E 

09, 

1 . 6  705000E 

05, 

4  2.95  760CCE 

C  3  , 

1 .68C7000E 

C  3 , 

8  •  8C36000E 

04, 

1.3581000b 

05, 

5-4  .  ?  7 !  ICOCfc 

C4. 

8 .223 JCOOE 

C4 , 

-1  • 9130C00E 

C9 , 

2. 6895000E 

C5, 

6-4.5C7400CE 

04, 

-2.97 19000E 

04, 

7.9432999E 

04, 

3. 1404000E 

05 » 

7-7.C111000E 

04, 

-5 .3676000E 

03. 

-1 • 9 1 JOCOOE 

C9  , 

4.9628C0UE 

C  5 , 

8- 3 . 525600  0  E 

C  5 , 

-8.972C999E 

04, 

-9.8 12  8999E 

04, 

6.8365CC0E 

04, 

9  3 . C34600CE 

C4 , 

3  .65  660  C  CE 

C4  , 

0. 

, 

9.  1  752999E 

02, 

1  3.61 7400CE 

C  8 , 

3.22550CCE 

07, 

3.5402000E 

C8  , 

1.658900CE 

C8/ 

DAT.MCP 

(1,1 

,1),I  =  153 

,  228)/ 

1-1.15C1000E 

C  6 , 

-2.321 100CE 

C6 , 

0. 

, 

2. 12  780002 

C  7, 

2  3.685400CE 

C  8 , 

3.246700CE 

C  7 , 

3.512 1000E 

C8 , 

1 . 6640000E 

08, 

3-3.E614C0CE 

C  6  , 

-3.2845C0CE 

C6  , 

C. 

, 

3.97C50GOE 

C  7 , 

4  3.330700CE 

ce, 

3.05040CCE 

07, 

3 • 7195000E 

08, 

1.6732000E 

08, 

5-1. 922700CE 

C  7 , 

2.924800GE 

C  7 , 

0. 

, 

1  •  188400CL 

08, 

6  6.1576999E 

08, 

3 . 0298CC0E 

C  7, 

3 .9904000  E 

C8 , 

2 • 555900CE 

ce, 

7-2 • 849000CE 

C  7 , 

4 .808 1999E 

07, 

C. 

, 

2.C49300GE 

08, 

8  1.57260CCE 

C8 » 

2.5620000E 

06, 

3 . S288000E 

C8, 

4.3566999E 

06, 

9-4.8845CC0E 

C  7 , 

-3 . 26560CCE 

C  7 , 

0. 

, 

4.11 540GCE 

C8 , 

1-2 .93  7600CE 

Cfc, 

-6.48290CCE 

07, 

2.3832000E 

C8, 

8. 7229999E 

C7 , 

2  1.25CC00CE 

03, 

1 . 80CC0C0E 

03, 

3 • COOCOOOE 

C4  r 

1.  10C0000E 

C4 , 

3  J.3CC0Q0CE 

03, 

6  i  TGCCOOOE 

02, 

2 • CCCCOOO  E 

C3 , 

5 .  CCCOOOOE 

C2 , 

4  4.500C00QE 

03, 

3 . 7000000E 

03, 

3.CCOCCOOE 

C4 , 

1.2250000E 

C4 , 

5  8.5CCOOOOE 

03, 

1 . 30CC000E 

C3, 

7.  CCCCCOO  E 

03, 

2.5000000E 

03, 

6  8.8CCC0C0E 

C3, 

8.50CCC0CE 

03, 

3 • OCOOOOOE 

04, 

1 . 3  500  000E 

C4, 

7  1 . fcCCOOOCE 

C4, 

3.C0C000CE 

03, 

9 • 5COCOOO  E 

03, 

1.4750000E 

C4 , 

8  1.22  5000CE 

C4  , 

1 .O5C0OCOE 

04, 

3.CCOOOOOE 

C4  , 

1.45C0000E 

C4, 

9  2.45C000CE 

C4 , 

1.5250C00E 

C4  , 

1 • 2C0C000E 

C4 , 

l.  77500GGE 

C4, 

1  1.60C0000E 

04, 

1.6CC00C0E 

04, 

3.000C000E 

04, 

1. 57500G0E 

04/ 

DATA(CP 

11,1 

,i),I  =  229 

,  240/ 

1  2. 70COOOCE 

C4  , 

2 • OOCOOCOE 

C4  , 

1.5COOOOOE 

04, 

2.6000000E 

04  , 

2  3.0COOOOCE 

C4  , 

1 . 80C0000E 

04, 

3.CCCOOOOE 

C4 , 

1.75C0000E 

C4, 

3  3 . CCCOOO CE 

C4 , 

3.00C000CE 

04, 

1.7500000E 

C4  , 

3. COCOOOOE 

C4  / 

OATA  NPtNSP 

,M,N 

,N5P1 ,hNNN/ 

4, 

8  ,  4 ,  5 , 

9 

,  6/ 

DATA  PSI 

/ 

1-0. 

, 

5 . OOCOOOOE 

-01, 

5 • CCCOOOO  E- 

01, 

5.  COCOOOOE- 

Cl, 

2-C. 

, 

5  .OOCOOCOE- 

-01, 

5.0CCOOOOE- 

Gi, 

0. 

3-C  • 

, 

0. 

-5. CCCOOOO E- 

Gi  , 

G. 

4-C. 

, 

0. 

, 

C. 

, 

0. 

5  C. 

, 

5.CCCC0CCE 

-01, 

5.CC00C00E- 

01, 

0. 

6  5.0CCG00CE 

-01, 

5.00C0C0CE- 

-01, 

5.CCC0GC0E- 

Cl , 

0. 

7  C. 

, 

0. 

-5.CCC0000E- 

Cl, 

0, 

8  C. 

, 

0. 

, 

C. 

, 

0. 

9  C. 

, 

0. 

- 1 • CCOOOOOE 

CO, 

0. 

1  C. 

, 

0. 

- 1 . CCOGQOQE 

CO, 

0.  ‘ 

2  G. 

, 

0. 

» 

1. COOCOOOE 

CO  , 

0. 

3  C. 

, 

0. 

-C. 

, 

c. 

4  0. 

, 

0. 

, 

C. 

, 

0. 

5  0  < 

, 

0. 

, 

C. 

, 

0. 

6  C. 

, 

0. 

-0. 

, 

0. 

7  C. 

, 

0. 

, 

c. 

, 

0. 

/ 

nn  onnfs  o  oooooooo 
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MMHIZL  FRtE  ENERGY 
NP  ARE  THE  NUMBER  OF  ELEMENTS 
NSP  ARE  THE  NUMBER  GF  SPECIES 
T  IS  ACCESS  FCR  KCPT=2,3 
KCPI=:  MFANS  P » T  GIVEN 

KCPT" 2  MEANS  P,F  GIVEN, PROP  IS  ENTHALPY 
KCP  T  =  3  MEANS  P,  S  GIVEN,  PROP  IS  ENTRCPY 
XII)  ARE  MCLE  FRACTIONS, YBAR=ICTAL  INITIAL  MCLES 
X ( NSP  1 1  =  I  • 

DC  23  1=1, NP 
A (NSP 1 ,  I ) =0  • 

DC  23  J=  I  ,NSP 

23  AINSPI,!  )=A(NSPi,  I)+A(J,I)*XlJ)*X(NSPl) 

RETURN 

ENTRY  GENR7X 
IF(KCPT.GT.OIGC  TL  26 
KCP T=-KOP  T 
CO  29  1=1, NP 
29  AINSPI,  I)=BEE(  I) 

26  WRlTE(6,llMxm,I  =  i,NP) 

26  CONTINUE 

DC  20  1=1, NP 
I F I  X ( I) .LE.O.IXII  )  =  .CQ1 
20  XLN(i)=ALCG(X(II) 

I  PH=  1 

2  FCRMATI3E1C.6, I  10) 

WRITE(6,13)T,P, PROP, KCPT 
13  FORMAT ( 3E20.B, I  10) 

INITIALIZE  lamua 

ZL AMUA-  1  . 

PLOG  =  ALOG(P  ) 

1999  GG  TC(20C7,2CCC,2CCC)  ,KGPT 
2CCC  IF(  IPF.EG.UGC  TO  2CC7 
LCFT=KCPT 

WILL  RETURN  IN  ETCH  THE  MIXTURE  ENTHALPY  FCR  KCPT 
1  =  2, ENTROPY  KCPT=3 
GO  TO  5 20C 

C216A  WRITE(6,11)ETCF,T 
2  16 A  CONTINUE 

IF IABS' (PRCP-ETCH)/PRCP) .L T . . COOCO i ) GO  TO  300 

FN2--F0LDH 

FOLCF=FNEWF 

FNEWH=PRCP-ETCH 

I F  C IPF.NE.2  )GC  TO  2CC2 

IN2  =  TCLC 

TCLL  =  T 

IF(FNEWH)2CC3,300,20CA 
2003  T  =  •>  95*  T 

GO  TO  2C07 
20 C A  T=1.C5*I 

GO  1C  2C07 

2002  IF( IPF.NE.3JGC  10  2CC8 

30  1C  FPKIMt=(FN£wF-FOLCF)/(T-TCLC) 

TN2=TCLC 
TCLC  =  T 

I=T-FNEWF/FPRIM£ 

GO  TO  200  7 
200  6  HN= I- TCLO 

nM=  I0LL-IN2 


1 

2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

1A 

15 

16 

17 

18 

19 

20 
21 
22 
23 
2A 

25 

26 

27 

28 
29 
3G 

31 

32 

33 
3A 

35 

36 

37 

38 

39 
AG 
A  1 
A2 
A3 
AA 
A5 
A6 
A  7 
At 
A  S 


ot 
5 : 

54 

55 
5< 
5 
51 
5‘ 
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i  I 

3 

?  | 

6  i 

7 

8  ! 


10 

11 

12 

13  I 


GN  =  FN/HM 

AM  =  9Ml  FNEWF-  ( l.HNI  *FCLCH*CN*FN2 ) 
bP=(2.*CMl.)*FNEWH-FCLCP*(  1  .*CN  )**2*0N*oN*FN2 
CP=FNEWF*(  l.  +  CN) 

CCN  =  8M*BM-4.*AW*Cf’> 

IF(CCN.LT.0.)GC  TC  3C1C 
CCMSCRT  (CCM 
0ENCM1  )  =  BHCCN 
DENCM2)  =  BP-CCN 

1F(ABS(CENCMI 1 ) )- AES ( CENCH i 2  ) ) ) 3CC2 , 3003 , 3004 
3CC2  BCT ICK=D£NCP ( 2 ) 

GC  TC  3005 
3C04  eCTTCP  =  OENCM  1) 

3C0  5  1N2=TCLC 
TCLC=T 

T=T-hN*2.*CM/BCTTCP 
2CC7  LCP  T=  1 

GC  TC  5407 
2163  KKK=  3 

GC  TC  4200 
4101  CCNTINCE 

C4  1 C  1  WRITE  (6, 11  )ElCF,Ffst*F,FPRlfcE,T  ,FCLCH,  ICLC,X 
2CC1  DC  12  1=1, NSP 
12  XL AST ( I )  =XLN ( I ) 

40  ITEH=C 
KCOM  =  C 

KKK  =  1 

64  KCtifS T 1=0 
LLL  =  3 
GC  TC  600 

210  DC  211  1=1, N 

211  FCLC(I)=f{II 

61  GC  7C(4C1/,4C2C),KKK 
4 1C2  LLL=2 

GC  TC  6 CO 

C45CC  HRITE(6,11)ZLAPCA,FQLC 
4500  CONTINUE 

I F  ( F  (  N  )  - FCL  C  ( N  ) )  7 C »  7C  »  7 1 

71  IF(F(N)  •LT..CCOCODGC  TC  65 
2LAPCA=ZLAMCA/2. 

KCLNT 1 =C 
KCUNT  =  KCUN  T ♦ 1 
KKK  =  2 

DC  75  1=1, NP 
FUI  =  FCLO(  l  ) 

75  XLM  I )  =XLAST (  1 ) 

F(N)  =  FCLD(M 
IF1KCLNT- 10)61,61, 72 

72  CALL  EKRCF(Gxl,KRLP, 131 
74  GC  TC  305 

70  KKK  =  1 

61  KCCNT1=KCCNT1+1 

1F(KC0NT1-2)E5,9Q,90 
90  2LAPDA=2.*2LAPDA 

IF(ZLAPCA.GT.1.)ZLAMCA=1. 

KCLM  1  =  0 

05  OC  66  1=1, NSP 

IFIABSI  (  XLM  I ) -XL  A  ST  i  t)  )  /  XLM  I)  )-. 00C1)B6, 86, 91 
66  CLNTK.OE 


tO 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 
7B 

79 

80 
61 
82 
83 
64 
85 
66 
87 
68 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
ICC 
101 
102 

103 

1 04 

105 

106 
1C  7 
ICB 

109 

1 10 
111 
112 

113 

114 

115 

1 16 
117 
l  ltt 


i 


IblHM-.LLLCCi  )SS,9S,<il 

51  IItR=lILR+l 

L  *RlTt(fc,ll)X 

L  :>AVb  LM  (  l  i  ,  1=1, NSP  ANC  flllihliN 

LL  t  S  1  =  1, NSP 
c  S  XL  A  b  T  l  1  )  =XLM  1  ) 

Ll  lfcS  1  =  1, A 
lcS  P ll C  (  1  1  =  F  I  1  ) 

IM  lTfcrt-JOcl  »S2,S2 

52  lALL  tKKlhUXl,KRlP,  12) 

SS  ClMlMt 

C  t>S  «K  1  I  b  (  fc  *  1  1  )  2LAPCA»X,XLN»XLAbT 
lHKCH.Et.nGC  TL  3LC 
l  PP  =  i  Ph*  1 

IF (  IPP-10) ISSS.iSSS,  ISSC 
l'jsa  CALL  PRftCMwXl  ,KKLM,  11  ) 

11  FlkPAT ( bc2C  .b  ) 

2 1 1  LlPT=KlFI+3 
GL  TL  52LC 
Jit  PHCP=bILP 
l  L  P  I  =  7 
Gl  TL  bSU7 
3 C  7  L  LP  I  = A 
3Cb  KCH  =  LCPT-3 
PR2=blCP 
Rb  1LRA 

Rl  1  7  LC  4CL1  J=i,l\P 
A  1  (  1 , J I =0  . 

LL  <«Cl2  1  =  1, P 
1  1  =  A  P  ♦  l 

SIC  2  Al(l,J)=Al(l,J)+X{ll)*etTA(i,J) 
sU  1  Alll,J)=X(J)+Ai(l,J) 

L  LF  FATK1X  IN  A1  ARRAY 

L  GbI  L l  1/YPAK  )/LLNX(L  )  l  A  IfcFPll) 

LL  SCCA  1  =  1  ,  A F 
TtPP(l)=G. 

CC  SCC3  J=N,i\SP 
J  J=  J-A  P 

SCI  2  IlRP(L)  =  IlRP(IMPSI{JJ,1,L)«X(J) 

SLOS  TbyP(L)  =  (T£PP(LH-xlL)*A(l,l))/A(ASPl,l) 

L  M*  PluP  UFUI/lLNXIL)  FCR  1  =  2,  RP  ANL  1=1, NP 

CC  SLC6  1  =  2  ,  A  P 
LL  S  0 1  fc  L  =  1  ,NP 
TbPPld  )=C. 

LL  SlLb  J=N,NSP 
J  J=  J--RP 

Sclb  UPPl(l)  =  JbFPl(lMFbl(JJ,l,l)*XU> 

SlCfc  A1 (  1 ,L )=A (L , I)*X( L )  +  TtPPl ( U-A(NiPi ,  1  )*TbPP( L) 
LLl  =  1 

S  C  l  7  tALL  ClbR /x ( A  1 ,LLNA , F  ,NP  ,N  ,  AT  ) 

SwL  LC  ACci  1  =  1, NP 

RC^l  XLMI)=ALMI)-21APlA*LINX(1> 

AclC  LL  ‘icC  1  1= 1 ,NP 
Sell  I  b  F  F  1  (  L  )  =  XLM  l  Mill) 

Cl  AC2b  1=R,NSP 
XLN (  1 ) =L  . 

J  J=  1-RP 

LL  Soil  L  =  1  ,  f\  P 

A  (;  J  1  XLMl)  =  XLMl)*JE*Pi(L)*titTA(JJ,L) 
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L  I"? 

4025 

XLM  1  )  =  XL.\(  I  >-UI  ) 

176 

1  f.\j 

DC  4C22  1  =  1, NSP 

179 

Lc  1 

IF ( XLM I ) )4C22,4C22,4C23 

16C 

l22 

C 

OCR  k  tC  T  ANY  MCLt  FKACIILNS  LARGER  IHAN  1 

161 

L  c  j 

4C23 

XLM1  )=-.Cl 

It  2 

L  ^  s 

4022 

CLMIME 

163 

L  c  b 

0 

CCMPLTE  X(I1  ANC  YBAK 

164 

Ub 

XINSP1  )  =  0 . 

165 

U1 

CL  405C  J  =  1 , N SP 

1 16 

1  <c  c 

X|J)  =  txPULMJl  ) 

16  7 

i  ib 

IF  (Ae  SI  XLM  J)  ).GT  .E8.  IX  (  J)=G. 

166 

i  jC 

4050 

X(NSPlJ=X(NSPI)*A(J»l)*vjJ) 

its 

1 3  1  1 

X(NSP1)=A(NSP1,U/X(NSP11 

ISC 

1  32 

GC  TC(4102,41C2,41C1),KNK 

191 

Hi 

0 

COMPUTES  F  FCR  NEfcTUN  RAFHSON 

192 

1  34  j 

c 

ALSL  STORES  LAPGESI  F  IN  F(Ni 

193 

lib 

oCC 

C0M=G • 

194 

1  it 

CC  6C1  1=1, NSP 

195 

1  37 

j  601 

L0M=C0M+X ( I ) 

1 9fc 

1  36 

F (  1  )=CUM-  1 . 

197 

i  34 

CG  6C3  1=2, NF 

196 

140  j 

C0M=C. 

199 

141 

CC  6C2  J  =  1 , N  S P 

2CC 

i.  42 

6C2 

DtM=DUM*A( J ,I)*X< J) 

2C1 

14  i  . 

603 

FI I)=CLM-A(NSP1,I  I/X4NSP11 

2C2 

144  1 

C 

F1NC  ThE  LARGEST  Fill 

203 

145 

F(M  =  ABS(F<  1)) 

2C4 

14c 

CC  6C5  1  =  2, NP 

205 

147 

E=A6S ( F ( I ) ) 

20e 

,  l 

146 

IF(fc-r<.\)  )  6C5  ,6C4 ,604 

2C7 

Il49 

604 

F  <  N  )  =  E 

2C6 

150 

6C5 

CCMINUE 

209 

151 

C 

kR  I TE 1 6 , 1 1 ) F 

210 

l  52 

GC  TC(4007,45CC,21C) ,LLL 

211 

1  b  3 

c 

FIRST  SUBSCRIPT  DEFINES  SPECIE  NUM3ER 

212 

154 

c 

SECOND  SUBSCRIPT  IS  FCR  A  RANGE  IN  TEMPERATURE 

213 

155 

c 

IFIRC  SUBSCRIPT  1=A 

214 

1  66 

c 

THIRC  SUBSCRIPT  2=B 

215 

15/ 

c 

TF IRC  SUBSCRIPT  3  =  C 

216 

1  5e 

c 

THIRC  SUBSCRIPT  4=C 

217 

154 

c 

THIRD  SUBSCRIPT  5  UPPER  I  IN  TEMP.  RANGE 

216 

1 1  0 

c 

S  =  A*LNT  +  BM*C 

219 

it  1 

c 

H=A*T*.5*B*T*T*C 

22C 

it2 

c 

CCNVtRT  MCLE  TC  MASS  FRACTIGNS 

221 

It  3 

c 

L  ARRAY  FCR  MASS  FRACTIONS 

222 

1  C  4 

520C 

EMM=C. 

223 

it5 

CC  bill  1=1, NSP 

224 

1  to 

5201 

EMM=EMM*X< I )*EM( I ) 

225 

16  7 

CC  52C2  1  =  1, NSP 

226 

1  fc  c 

5202 

/  (I1  =  X( I ) *  £ “ I  I ) /EMM 

111 

its 

C 

WKITE(6,11)2*EM,EMM 

lib 

l  7  0 

540  7 

T LLG  =  ALLG l T  ) 

229 

1  71 

E  TCH  =  C . 

230 

1  7  <L 

CC  54C1  1=1, NSP 

231 

1  7  j 

DC  5hC2  j J= 1 ,NNNN 

232 

l  74 

IF (CPI  I, JJ,5)-T ) 5402, 5403, 54C3 

233 

1  lb 

5402 

CCM  INUE 

234 

1  7  o 

CALL  ERRCR( T ,  1 , 10 i 

235 

1  7  / 

5403 

GC  IC I 5406, 5405, 54C4 ,5405, 54C4.54C5, 5404) ,LGPT 

236 

5405  hi « (CPU*  JJ,  1)*.5*CP4I.JJ,2)*T)*T*CP4  l,  JJ,4) 

C  MRITE46, 113M»  JJ,HI 

ETCF=ETCh*Z  4  I  ) *h I 
GC  TC  5401 

5404  SI=CP(I,JJ,1)*TL0G+CPI if JJ,2‘*T+CP4I ,JJ,3) 

C  )«RiTE  16*  1131  I*JJ«  SI 

ETCF  =  ETCMX4  I  )* IS  I/M  I  I l-XlNiEI-PLCGI /EMM 
GC  TC  5401 

540  fc  CUi=CP(  I  tJJ,ll-CP4I  ,JJ,3I-.5*CPI  If  JJ,2)  *T+CPI  It  JJt4)/ 
1T-CPII,JJ,1)*TLCG 
C(I)*Cii)/M4I)+PLCG 
U  MHITE46, 1131  it JJtCI I  I 

5401  CONTINUE 
C  IM  FORMAT (21 IC.3E2C. El 

GC  TGI  2  16 3, 2  164, 2 164, 206, 30  5,305,3071 tCQPT 
ENC 

SI  SCLEM7 

SUBROUTINE  CLEM7X4B,X,C,M,M1,AT) 

DIMENSION  BIM,M),X(M),DIM1I,AT4M,MI) 

C  EOLATIONS  ARE  OF  Th£  FORM  BX  =  D 

DC  17  1  =  1, M 
17  XIII=C.G 

CO  2CC  1=1, M 

200  A T (  I  ,M1 )  =  C( I ) 

CO  20  1  1  =  1, M 
DC  2C1  J= 1 , M 

20 1  ATI  I  ,J)  =  B( I ,J) 

DO  32  N=l,M 
C*AT4N,N) 

I  l  =  C 

’  DC  $  I =N , M 

IF  (A6SIAT  (I,N))-AES(C))9,S,8 
E  C  =  A 1 1 1  ,N ) 

II  =  1 

S  CONTINUE 

IF!  I T-N  )  7 , 7 , 70 

70  DO  71  J=N,M1 
TEMP=AT I N,J) 

ATIN, J)=AT< I T , J ) 

71  ATI  IT , J)  =  TEMP 
7  CO  10  1  =  1, Ml 

1C  ATIN,I)=ATIN,IJ/0 
IFIM-N)50,5C, 1€ 

IE  N1  =  N+1 

CC  30  I  =  N  1 , M 
C  =  A  T I  I , M 
CC  20  J  =  N  ,M1 

30  AT(I,J)=AI4I,J|-ATIN,JI*C 
32  CCNTINUE 
50  X  I H  1  = A 1 1 M ,M  + 1 J 
DC  65  N=2,M 
NR=M*1-N 
C  =  A T I  NR , M*  1 ) 

DC  60  I  =  NR  ,M 
60  C=G-AT (NR  ,  I )*X( I) 

65  XINR)=C/AT(NR,NRI 
RETURN 
ENC 


412  CARCS  PRINTEC 


237 

238 

239 

240 

241 

242 

243 

244 
243 

246 

247 
246 
249 
25C 
231 

252 

253 

254 

255 

256 

257 

258 

259 
26C 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
261 
282 

283 

284 

285 
266 

287 

288 

289 

290 

291 

292 

293 

294 


APPENDIX  C 

LISTING  OF  SUGGESTED 
SUBROUTINE  ERROR 


C-l 


iiehTC  LIST  LURCH* 

SCtiKLCT  I  N  t  fcRRlRU*K»Ii 
WHlTtlOtiOI 

IC  F CRH  AT  (6H.ERRCH»  f  3) 

J=  I-S 

oC  TCU,4»J»4J,J 
1  y»K  I  T  b  (  fc  f  i  1 )  C  » K 

11  f-LKMMI  l<*h(JTtHFcKAIUkt  =  1 1  i  .  i  t  IS*-  F  I  T  EUR  bPLClL  NO.  htVhMf 
1  ) 

GC  TC  7 
l  wKlTt(6*lz!j 

VI  PukRA  I  (431-OICTcR  twUlLlERlCR  ileKATICNS  tXCfcfcC  1G 
Gc  IL  7 
i  nKllt(bfli) 

Li  FtRHAT ( A^bCPREaSoKt  TtRPEkATCKL  IlfckAlIUNh  tXCEEC  3c 
GL  IC  7 
4  itRITcUf  14) 

14  PCkR  A  T  1  l'jhCLARCA  CCT  1G  TIRES 
7  RcTCHk 
ENL 


VAL  II 


) 

) 

) 


CAHGS  PRIMEC 


APPENDIX  D 

SAMPLE  INPUT  FOR  THE 
CHEMISTRY  EQUILIBRIUM  GENERATOR  PROGRAM 


c 


WORD  2 


3  1 

0 

D 

0 

D 

WORD  5 

WORD  6 

WORD  T 

WORD  8 

D 

D 

D 

D 

D  D  a  0  0  D  0 

WORD  I  WORD  2  WORD  3  WORD  4  WORD  5  WORD  6  WORD 

0  0  Q  D  D  D 


D  D 

WORD  8 

D  U 


TWFnyoTTo  OTOl  o:o  o  c  ojc  o;o  o  o  sio  o  o  s  o  Fo  o  6  q]o  o  o  c!o  oio  o  o  o|o  i  c(o  oio » o  o I I  o  o|o  r:n  o  o  o]5  dT L , i  Fa' 0  0  ofa  c  c  0 

1  )| )  <  >  >i  1  1  1  10  11  tail  H  is  ujr  ii  :i  x  :>  :j,r.  ji ss  3  :i  n  n  »  !'  uln  «  is  *  1 3*  :• «  •'  u  is  *l<i  a  is  ssisi  s.)1)  si  ;s  «jsi  ss  ss  n  >i  iAj  w  is  fill  t»  m  .0  n  ss|’i  ’  1*  i|  t 's  * 

1  ],;  1 1 1111 1 1)  n  1 1 1.1 11 1 1  i>i  1 1  i!i  1 1 1 1  in  1  in  1 1 1 1  i|i  1 !  in  1  1  ih  ;‘i  1 1 1.1 1 1 1 1  i|i  1 1  i|i  1 1 1  1  i,i  i ; ;  1 1 1  i 

2  2|2  2  2  2|2  2  2  2  2  2j?  2  2  2|?  2  7  ?  2  ;j?  2  2  2*2  7  2  2  2  2(2  2  2  2|2  2  2  2  2  2I2  2  2  2*2  2  2  2  7  2|2  2  2  2>2  2  2  2  2  2l2  2  2  2*2  2  7  2  2  2j2  2  2  2j2  7  7  ? 

3  3^3  3  3  3|3  3  3  3  3  3|3  3  3  3|3  3  3  3  3  3|3  3  3  3^3  3  3  3  J  3J3  3  3  3|3  3  3  3  J  3*3  3  3  jjl  3  3  3  3  J|3  3  3  3*3  3  3  3  3  3*3  3  3  3*3  3  3  3,3  3*3  3  3  3>3  3  3  3 

4  4^  4  4  4|4  4  4  4  4  4j4  4  4  411  4  4  4  4  4|4  4  4  4  4  4  4  4  4  4^  4  4  4|4  4  4  4  4  4*4  4  4  4|4  4  4  4  4  4j4  I  4  1*4  »  «  4  <  tja  4  4  4^  4  4  4  4  1*4  4  4  lla  4  4  4 

5  J|5  5  5  S*S  5  5  S  5  5jS  S  5  5|s  5  5  5  5  5I5  5  5  5|5  5  5  5  I  SlS  5  5  S|5  5  5  5  5  5  5  5  5  5|5  5  5  5  5  5|5  5  5  5|5  5  5  5  5  5^5  5  5  5jb  i  b  5  5  5,5  5  b  5,5  5  5  b. 

I5  5*5  6  5  eJs  8  6  8  6  t|i  5  S  (It  6  6  6  6  f*3  I  (  t|(  S  (  (  1 1*4  4  6  fU  t  8  6  5  6*6  6  6  (|l  I  (  (  6  6*1  6  (  (jfi  6  «  5  Cp  t  i  S'i  B  C  5  6  t|8  5  0  «>«  S  5  * 

L  jjn  1  ip  ?  1 1 1  ?l?  j ;  >*3  7 ;  ?ii  »*n » ;|i  1 1 1  n*n  i  ?*?  1 1 1 1  ?!? ;  ?  il?  1 1 1 1  ?*»  n  il?  1 1 1 1  ?!i  nAiiinlinilim 

a  ala  a  a  ala  a  a  a  a  a*a  a  a  a*a  a  a  i|a  a|»  a  a  ala  a  a  a  a  a  a  a  a  a*a  a  a  a  a  a{a  a|a*a  a  a  a  a  a*aa  a  a'a  a  a  a  a  a*a  a  a  a!»  a  a  a  a  a*a  a  a  a*a  a  a  t 

|  I  !  *  *  I  I  *  1  *  II  !  *  I  I 

9  9  9  11  9  9  9  9  9  9  "  9,9  9  0  9i9  9  9  3  3  9,S  !  9  9  9  9  8  9  9  9i9  9  8  9|9  9  !  9  9  9]9  !  9  9i9  9  9  9  9  9.9  9  9  9(9  9  9  9  9  9  3  9  9  9.9  9  9  9  9  9(!  M  9|9 1  9  9, 

1  2!  j  ;  5  il  1  t  9  !«ti  »»»’j  :m  a(lij  :•  ;3rt*7!  »  n  mvnnnv  Hjj  h  n  *  v  *  a  «n  «?o  •««!  <r«?  «  wsdji  aHu  »»  als? » *idti  dnMMuUrutt  win  irn  .*  .Vi;  »ui 


TAQCO  r  -.6424 


D  0  0 

WORD  I  I  WORD  2  WORD  3  WORD  A  I  WORD  5  I  WORD  6 

D  0  0  0  0 


WORD  7  WORD  8 

0  0 


il'O 

0  B  i  8  0  0 

0  1  Hi 

0  0  o' 

1  ill 

<  s  ih  •  1 

tolti  inn  14 15  it 

1  'l' 

1 1  ip  1 1 1 1  ri 

1 1 1 

nn 

2  2  2|2  2  2 

2  2  2,2 

2  2  2 

3  3I3 

3  3  3(3  3  3 

»v 

3  3  3 

4  4  4 

4  4  4|4  4  4 

4  44 

5  5|5 

5  5  sU  5  5 

5  5  5,5 

6  ill 

1 

1 1  l|i  6  1 

liili 

i  i  i 

1 

7  7  7 

7  7  7(  7  7  7 

7  7  7b 

1 

1  1  1 

1  iii 

1 1  ilm 

10  Ji 

1  *1 

S  9  9 

1  iii 

S  1 1 1 1 S 

4  S  III  1  1 

I  I  8.S  S  9  Si 

1«U  12113  14  15  III 

g  0(0  v  0  0|S  0  0  0  0  o'l  0  0  0|M  0  D  0  0|9  ■  D  l|B  ■  D  U  B  BID  I  DIO  0  e  0  0  0(0  B  0  0(0  S  B  0  0  0*0  S  0  gi 

ii  HJ3  ii  is  it'ii  it  a  a  n  »|a  h  is  y  a  a  « ti  tAi  «  is  *'ti  a  ti  a  si  situ  st » a|si  a  so  u  it  h'o  «  a  *\>  a  u  it  it  n|’i  it  is  hi 

1  v\  1 1  ill  1 1 1  i-i|i  1 1  ih  1 1 1 1  ill  1 1  i|i  1 1 1 1 1 1 1 1  i|i  1 1 1 1  i|i  1 1  iii  1 1 1 1  i|i  1 1 ; 

2  2  2  2  2  2^2  2  2  2  2  2I2  2  2  2.  2  2  2  2  2  2^  2  2  2^  2  2  2  2  2i2  2  2  2|2  2  2  2  2  2I2  2  2  2^2  2  2  2  2  2|2  2  2  2| 

I  .1  .1  I  .1  ' 


1  in  1  mm  1 1  -i| 

!  2  2  2  2  2^2  2  2  2  2  2I 

1  I  I 

nil  n  3 1 1 1 1 1 11 

1  I  I 

>  4)4  4  4  4  4  4  4  4  4  4. 
i  S|S  5  5  5|S  5  5  5  S  Si 

1  sis  I  I  III  I  I  I  I  |l 

I  1 

1  n  1  j  ;|m  n  ? 

1 1 1 1 1  il|  1 1 1 1 1 

I  I  I 

1  s.im.n  1  m, 

1  »n  h  is  Art  21  i*  in  tdr 


1 3  3|3  3  3  3  3  3>3  3  3  3,3  3  3 : 

*  I 

14  4(4  4  4  4  4  14  4  4  4(4  4  4  1 

1 5  S|S  S  5  S  5  Js  5  5  S|5  5  5 ! 

1 1  sic  6  I  fi  S  i|s  S  S  sis  S  1 1 

iJiiiiiiniiliih 
I  1  I 

iiijimiijiiiijmi 

ii  si  1  n  1  nm  i.mi 

1  kWh  a  ado  du  ««*<!««> 


oTTI Him  0(0  e  s  0  0  o'o  a  0  0)0  0  o  n 

si  a  so  u  si  a'o  m  is «  v  a  u  it  11  »|’i  it  is  nin  11  s » 

1  1  1  1  1  1I1  1  1  iii  1  1  1  1  iii  1  t  :  1  1  1 

2  2  2  2  2I2  2  2  2b  2  2  2  2  2|2  2  2  2i2  2  2 

|  ’ 

3  3  3  3  3  3  3  3  3  3  3  3  3  3  ]l]  3  3  ]|j  3  3 

I  I  ll 

4  4  4  4  4  4  4  4  4.4  4  4  4  4  4'4  4  4  4>4  4  4 

I  1  II 

5  5  5  5  5,5  5  5  S|S  S  S  5  5  5|5  S  5  5|5  S  5 

I  I  i  6  ill  i  6  lls  I  6  6  i  l|i  i  i  i|l  i  i 

m  n\m  ri  1 1 1  n\i  m\n  1 

i  1  i  i  ill  1 1 1 1 1 1 1 1  i'i  1 1  JiOi 


?|i  in'11 


TiBCO  *6642 


9 1,9  9 1  9,i  9  9  9  9  Si  S  S  1,1  1 1  19  19  S  1 1,9  991 

II  du  W  »  Mis?  X  99  Mil  UU  MB  MW  M  IS  70|*l  )|7|  74  ?J  Ml/  71 7*  N 


WORD  e 


HBMUWBS 


i rm i.i ii  i 


3  3<3  3  3  3(3  3  3  3 


2  2  2  2  2  2,2  2  2  2 


l[lHriI|[IIfl[IUBII|fTT|ll 


si  ew  M  »  » v  m  m  Tipi  7rjn  73  n  it  j»7  ti  ti  id 

1  ill  i  1  ill  i  i  i  i  1,1  1  i  :Ji  i  i  i 

2  2I2  2  2  2^  2  2  2  2  2|2  2  2  2|2  2  2  2 

1  1*3  3  3  1^3  3  3  ?  3  3l3  3  3  ill  2 1 1 

4  4  4  4  4  4,4  4  4  4  4  4*4  4  4  4*4  4  4  4 

I  1  II 

|  5,5  S  5  S|S  5  5  S  5  5,5  5  5  5.S  S  S  S 


1  7|7  7  7  I'l  1  1  l\l 


iii  ill 


D-3 


CARD 

NO. 


STANOARD  8  WORD  CARD 


